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Abstract 


The dynamics of excited electrons and holes in single layer (SL) M0S2 have so far been 
difficult to disentangle from the excitons that dominate the optical response of this material. 
Here, we use time- and angle-resolved photoemission spectroscopy for a SL of M0S2 on a 
metallic substrate to directly measure the excited free carriers. This allows us to ascertain a 
direct quasipaiticle band gap of 1.95 eV and determine an ultrafast (50 fs) extraction of excited 
free caniers via the metal in contact with the SL M0S2. This process is of key importance for 
optoelectronic applications that rely on separated free carriers rather than excitons. 

KEYWORDS: transition metal dichalcogenides, M0S2, free carriers, excitons, time- and angle- 
resolved photoemission spectroscopy. 
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The optical properties of semiconductors are completely specified by the dynamics of free 
carriers and bound excitons, both stemming from electron-hole pair excitations.- In bulk semi¬ 
conductors such as GaAs- and two-dimensional ( 2 D) quantum well systems,-^ these excitations 
have been studied for over half a century. Recently, genuine 2 D semiconductors such as single 
layer (SL) M0S2 have been isolated.-^ In these materials the low dimensionality contributes to the 
presence of strongly bound excitons- with remarkable properties. 

The current understanding of the optical properties and excited carrier dynamics in SL and 
few layer M0S2 relies largely on a series of photoluminescence and differential absorption ex¬ 
periments. The electron-hole pair dynamics has been investigated through the measured exciton 
lines, which are the most prominent features in such experiments.— This approach has been re¬ 
markably successful leading to, for example, the demonstration of a valley polarization between 
optically induced excitons associated with the degenerate conduction band (CB) valleys in the 
material’s band structure.— Time-resolved studies have revealed valuable insights into ex¬ 
citon lifetimes, including relaxation dynamics that evolve over femto- to picosecond timescales 
due to scattering processes involving lattice vibrations-^ and defects.— However, an important 
technological challenge remains: to separate the bound electron-hole pairs into free carriers that 
can then be extracted via metallic electrodes. This approach has recently been explored in com¬ 
plex 2 D material heterostructures consisting of graphene and SL transition metal dichalcogenides 
(TMDCs),— paving the way towards efficient photovoltaic devices with improved absorption 
characteristics 

While the all-optical studies of SL M0S2 have unveiled the exciton physics in this material, 
the strong exciton-related features in the spectra mask the free carrier contribution to the optical 
properties almost entirely.— We overcome this issue by assessing the free carrier dynamics in SL 
M0S2 with time- and angle-resolved photoemission spectroscopy (TR-ARPES). For this study, 
we employ a model system consisting of high quality SL M0S2 epitaxially grown on a Au(lll) 
substrate (see sketch in Fig. [Ua)).— Not only does this system allow us to quantify the direct band 
gap of free carriers in a SF of M0S2 on Au(l 11 ), it permits us to measure free carrier lifetimes in 
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SL M 0 S 2 near a metal contact - a configuration that is relevant to real device architectures.-^^jM 

Using TR-ARPES to unravel the carrier dynamics near arbitrary points of a material’s Brillouin 
zone (BZ) has only recently become possible with the advent of high photon energy laser light 
via high harmonic generation.—-— This is especially relevant for electronic states near the BZ 
boundary, such as the Dirac cone in graphene^>^ or, indeed, the CB and valence band (VB) 
extrema of SL M 0 S 2 , which are sketched in Fig. [IIb)-[Ttd). In the TR-ARPES experiment, the 
sample is first excited by an ultrashort tunable pump pulse, such that optical pumping can be 
achieved with energies below (Fig. [Hb)), at (Fig. [He)) and above (Fig. [ltd)) the quasiparticle 
band gap. The system is then probed along the F — K direction in the BZ while in the excited state 
by photoemission using a time-delayed 25 eV probe pulse, as sketched in Fig. Uta). Both pump 
and probe pulses have a duration of 30 fs. 

The electronic structure of our sample at equilibrium is characterized by sharp dispersing VB 
features of SL M 0 S 2 with spin-orbit split bands around the VB maximum (VBM) at K, as seen in 
the ARPES data in Fig. (He).— The location of the projected bulk band continuum and bulk gaps 
of the underlying Au(l 11) have been added, because these will play a crucial role for the dynamics 
observed in the three different optical pumping regimes. 

The non-equilibrium carrier distribution is displayed in Fig. [Ilf)-[Ilh), which show the intensity 
difference between the spectrum immediately after the optical excitation (30 fs) and a spectrum 
taken before the arrival of the pump pulse. In these difference plots, red and blue indicate an 
increase and decrease of spectral intensity, respectively. As expected, such a difference spectrum 
is rather featureless when the system is pumped with 1.6 eV, an energy well below the band gap 
(Fig- iHf))- When the excitation energy is increased to 2.0 eV, we observe a clear population of 
the CB states at K, as seen in Fig. Htg). This indicates that our pump energy exceeds the band 
gap and excites free carriers into the CB, a point we return to below. This population strongly 
decreases when we excite with 3.0 eV photons (Fig. H^h)). For this excitation, one might expect 
a population of electronic states at higher energies in the CB and an eventual cascading down 
to the CB minimum (CBM). In this scenario one would either expect that the intensity at the 
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Figure 1 : Tuning pump-probe ARPES experiments near the direet band gap of SL M0S2: (a) 
Sehematie of our photoemission experiment. Eleetrons (blaek spheres) exeited by a pump pulse 
tuned aeross the visible to infrared range are photoemitted from epitaxial SE M0S2 on Au(lll) 
using a 25 eV probe pulse and deteeted along the high symmetry direetion E — K in the hexagonal 
BZ. (b)-(d) Diagrams of eleetron-hole pair exeitations around K for the studied regimes of optieal 
exeitation: (b) below the CB onset, (e) at the band gap energy and (d) above the CB onset. Holes 
are sketehed as orange spheres. The pump pulse energies and eorresponding wavelengths used in 
the experiments are stated, (e) Oeeupied band strueture at equilibrium measured by synehrotron 
radiation, (f)-(h) Intensity differenee obtained by subtraeting the speetra at the peak of the optieal 
exeitation from a speetrum obtained before the arrival of the pump pulse. The orange lines in (e) 
and full blaek lines in (f)-(h) outline the boundaries between the projeeted bulk band eontinuum 
and gaps of Au(lll). The dashed blaek eurves in (f)-(h) are the fitted dispersion from (e). The 
optieal exeitation energies in (f)-(h) eorrespond to those given direetly above in (b)-(d). 


CBM would inerease at a later point in time or that the rate of deeay of eleetrons in the CBM 
would slow down. We are only able to observe weak signals at the CBM with the 3.0 eV pump 
pulse, whieh we attribute to the dual eonstraints of phase spaee restrietions imposed by energy and 
momentum eonservation rules for photo exeitation of M0S2 at this energy and a shorter lifetime 
of the M0S2 states exeited within the eontinuum of metal bulk states. The latter is fully eonsistent 
with the inerease of linewidth, and thus deerease of lifetime, of the VB states as they eross into the 
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Figure 2: Snapshots of excited free carriers around the K point: (a)-(e) Intensity difference with re¬ 
spect to a spectrum obtained before optical excitation {t < 0) at the time delays marked on the time 
line. Red contrast can approximately be interpreted as excited electrons while blue corresponds to 
excited holes. The position of the band edges and Fermi level are shown by dashed horizontal lines 
with the energy differences given in electron volts in (e). (f) Energy distribution curves (EDCs) of 
the intensity difference (grey curves) in the VB and CB regions at a time delay of 30 fs. The EDCs 
have been binned over a range of ±0.01 with respect to the given kx values. Black curves are 
fits to Gaussian functions. Blue and red parabolic traces follow the peak maxima in the VB and 
CB, respectively. The uncertainties on the band offsets in (e) and the stated band gap value {Eq) in 
(f) are ±0.05 eV. All data were taken at an excitation energy of 2.0 eV. 


continuum of metal bulk states, which can be seen in the equilibrium ARPES data in Eig. (He) and 
in the more detailed analysis of these states in Ref.— 

The data in Eig. [TJg) establish that optical pumping with 2.0 eV photons is required to resolve 
the CBM around K and thus to observe the free carrier dynamics of SE M 0 S 2 . We have therefore 
used this excitation energy while recording the time-resolved series of snapshots of this region of 
the BZ. Snapshots of this series are shown in Eig. [2] while the quantitative intensity development 
is discussed later in connection with Eig. HI Excited electrons and holes are already distributed 
around the VB and CB edges during the optical excitation phase at t = 0, i.e. at the centre of the 
pump pulse (see Eig. I^a)). This is ascribed to rapid carrier-carrier scattering. The peak excitation 
is reached at t = 30 fs in Eig. I^b) (see also Eig. H]), where the full parabolic dispersions of the 
band extrema emerge due to the increased population of carriers. This then rapidly relaxes with 
the overall intensity difference already halved at t = 60 fs in Eig. [2tc) and approaching the noise 
level at t = 90 fs in Eig. [2td). In Eig. [2te), at t = 120 fs, no excited carriers can be seen. 
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The parabolic dispersion around the VB and CB edges can be extracted from the peak intensity 
difference at t = 30 fs using the analysis introduced in Fig. [2f). This allows us to extract values for 
the VB and CB offsets with respect to the Fermi level, as indicated in Fig. [2te)- We can therefore 
directly measure a free carrier direct quasiparticle band gap of (1.95 ±0.05) eV, which is marked 
in Fig. [2tf). This value is different from a previous estimate by high resolution ARPES using the 
standard approach of populating the CB by alkali-doping.-2^>2^ We ascribe this to the possibility of a 
severe band distortion caused by alkali-MoS 2 bonding, inhomogeneous doping or band tailing 
effects. The observed gap is substantially smaller than a theoretically estimated value of 2.8 eV for 
free-standing M 0 S 2 .— This can be explained by a strong renormalization of the band gap due to 
the metallic substrate. Indeed, similar renormalization effects, although substantially weaker, have 
already been observed for 2D TMDCs grown on graphene on silicon carbide and graphite.—*^ 
Note that the quasiparticle gap that we determine for excited free carriers here is a fundamentally 
different quantity from the optical gap that derives from the energy of the A exciton line. The 
size of the optical gap is around 1.9 eV for exfoliated M 0 S 2 placed on insulating substrates such 
as Si02/Si— and hBN/Si02/Si.— In our case, the close proximity of SL M 0 S 2 to the underlying 
metallic substrate is expected to cause a significant screening of the excitons.— This reduces the 
exciton binding energy, which may explain why the measured size of the quasiparticle gap for 
M 0 S 2 on Au(l 11) nearly coincides with the optical gap for M 0 S 2 on Si02/Si. 

A more detailed analysis of the photoemission intensity for a pump excitation with 2.0 eV 
photons is presented in Fig. Oa). The black and orange curves are energy distribution curves 
(EDCs) of the photoemission intensity at t < 0 and t = 30 fs, respectively, integrated over the 
momentum range that contains the CB and VB edges. The EDC at t < 0 was obtained from a 
spectrum acquired 200 fs before the arrival of the optical excitation. Pumping the system causes the 
t = 30 fs spectrum to be slightly different from the equilibrium situation at t < 0, with a depletion 
of intensity in the VB part and a small rise of intensity above the noise level in the CB region. 
This enhanced population of electrons around the CB edge and the corresponding holes around 
the VB edge are clearly seen in the difference between these EDCs (grey curve in Eig. Oa)) in a 
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Figure 3: Analysis of the photoemission intensity for exeitations at the gap energy (2.0 eV) and 
above (3.0 eV): (a) EDCs at 2.0 eV optieal pumping, extraeted by integrating the absolute photoe¬ 
mission intensity over ±0.3 A ^ with respeet to K. An equilibrium EDC obtained before optieal 
exeitation {t < 0) and an EDC at the peak exeitation {t = 30 fs) are shown for eomparison. The 
CB (VB) edge is identified via a small surplus (depletion) of intensity, which can also be seen by 
a peak (dip) in the plotted difference between the two EDCs (grey curve). The excitation process 
around the VB and CB edges is defined in the upper right diagram, (b) Plots of EDCs integrated 
over ±0.3 A^^ with respect to K as in (a) but at a photon energy of 3.0 eV, well above the gap 
energy (see upper right diagram for definition of this process). A shift on the order of 10 meV 
between the two EDCs is highlighted by a double-headed arrow. This is clearly visible via an 
increase in the difference curve around a binding energy of 1.0 eV. 


similar way as in the corresponding Eig. [Ug). The same analysis for the 3.0 eV optical pumping 
regime is shown in Eig. Ob). The analysis reveals that the intensity difference signal in the VB 
region is instead dominated by a small shift of the EDC at t = 30 fs to lower binding energies, 
compared to the EDC taken before excitation. This yields a substantial increase (decrease) in the 
difference curve on the lower (higher) binding energy side of the VB peak, which is also visible 
in the entire M 0 S 2 VB in Eig. Oh). While the changes in the intensity difference seem dramatic, 
this shift is merely on the order of 10 meV as seen in the EDCs in Eig. Ob). We can rule out that 
the shift is caused by a surface photovoltage effect or by a photoelectron space charge above the 
sample surface because only the the M 0 S 2 states shift. Both aforementioned effects would have an 
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Figure 4: Efficient extraction of excited free carriers by a metal contact, (a)-(b) Time dependence 
of the summed intensity difference (filled dots) in (a) the CB and (b) VB regions of the data. Fits 
to functions consisting of a rising exponential function and a single exponential decay broadened 
by the time resolution are shown by black lines. The stated values are the time constants of the 
exponential decays, (c) Interpretation of the observed dynamics. Electrons (black circles) excited 
into the CB of SL M 0 S 2 recombine with holes (orange circles) in the Au(lll) bulk continuum 
(orange shaded region) by ultrafast Auger processes (black curled arrows). The same type of 
processes (orange curled arrows) are responsible for the recombination of holes in the VB of M 0 S 2 
with electrons in Au(l 11). 


impact on the photoelectron propagation in vacuum and thereby lead to shifts of both M 0 S 2 and 
Au states.—*^ 

For the data taken with a 2.0 eV pump pulse, the shift appears insignificant in the intensity 
difference, and we can reliably assume that the time dependence of the intensity in the CB and 
VB regions originates entirely from the electron-hole pair dynamics induced by the pump pulse. 
This is studied in further detail in Fig. |4ta)-l4tb), showing the time-dependent integrated CB and 
VB intensity, respectively. In both cases we find the decay part to be well-described by a single 
exponential function with a time constant of 50 fs, which is at the limit of the time-resolution of 
our experiment. The free carrier dynamics observed here are several orders of magnitude faster 
than the dynamics reported for bulk GaAs which is also a direct band gap semiconductor.— The 


9 








50 fs time-scale is too short for lattice excitations to take place. The observed ultrafast dynamics 
can therefore only be explained by Auger-type electron-electron interactions directly involving the 
metal contact, as sketched in Fig. I5c). Excited electrons in the CBM of SL M 0 S 2 can easily 
find a hole in the nearby bulk continuum of states in the metal while simultaneously exciting an 
electron-hole pair in the Fermi sea of the metal. The reverse of these processes accounts for the 
hole dynamics in the VB. 

Finally we note that the free carrier dynamics we have observed here are also much faster than 
any other free carrier decay mechanisms - for example, exciton condensation, recombination in 
trap states originating from defects or radiative recombination. This suggests that extremely effi¬ 
cient free carrier extraction is already present near gold electrodes in contact with the SF M 0 S 2 . 
This type of material combination thus comprises a route towards exploiting the optically gener¬ 
ated free carriers in novel 2D semiconductor systems. 


Methods. The SF epitaxial M 0 S 2 samples were grown by a physical vapor deposition (PVD) 
procedure consisting of Mo evaporation on a clean single crystal Au(l 11) substrate in an ultra-high 
vacuum (UHV) chamber with a low pressure of H 2 S gas, followed by annealing in the H 2 S atmo¬ 
sphere. The growth was cycled to provide large-area SF M 0 S 2 as described in further details in 
Ref.— High resolution equilibrium ARPES data were obtained by transferring the sample through 
air to the SGM-3 UHV end-station of the synchrotron radiation source ASTRID2 in Aarhus, Den¬ 
mark. The sample was annealed to 500 K to remove adsorbed species and to give rise to sharp 
M 0 S 2 electronic states in ARPES spectra (see Fig. [He)). The ARPES data were collected with the 
sample temperature kept at 70 K. A photon energy of 25 eV was selected in order to be consistent 
with the TR-ARPES data. The total energy- and angular-resolution amounted to 20 meV and 0.2 °, 
respectively. 

The TR-ARPES data were acquired at the Artemis facility, Rutherford Appleton Faboratory in 
Harwell, UK.— Samples were transported in an evacuated tube pumped down below 1 • 10~^ mbar. 
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Once the samples were placed in the TR-ARPES end-station they were annealed to 500 K to 
remove any adsorbed surface contaminants. The sample temperature was kept at 50 K using a 
liquid helium eryostat during measurements. A 1 kHz Ti:sapphire amplified laser system with a 
wavelength of 785 nm, a pulse duration of 30 fs and an energy per pulse of 12 mJ was used to 
generate pump and probe pulses. In order to measure the band strueture throughout the entire first 
BZ of SL M 0 S 2 , probe pulses of hv = 25 eV, corresponding to the 15th harmonie of the laser 
fundamental, were generated by foeusing a part of the laser energy on a pulsed jet of argon gas. 
The remaining laser energy was used to drive an optical parametric amplifier (HE-Topas) followed 
by a frequeney mixing stage, whieh enabled us to obtain tunable pump pulses with wavelengths 
eentered at 408 nm, 615 nm and 784 nm. The applied pump fiuenees were kept in the range 1.4 
to 4.8 mJ/em^, and were always optimized to avoid pump-induced spaee-eharge broadening of the 
acquired speetra.— We used a eonfiguration where the pump pulse was 5-polarized and the probe 
pulse was p-polarized. The time delays of pump and probe pulses were varied using a mechanical 
delay line. In order to aehieve a satisfaetory signal to noise ratio in the pump-probe spectra total 
acquisition times on the order of 17 hours per dataset were neeessary. The energy-, angular- and 
time-resolution were 400 meV, 0.3° and 40 fs, respeetively. The intensity differenee in Eig. |4ta) 
was summed over a binding energy range from —0.95 eV to —0.20 eV and a momentum range 
from 1.0 A ^ to 1.5 A“^. In Eig. |5b) the sum was performed over the same momentum range but 
over a binding energy range from 1.00 eV to 1.75 eV. The exponential funetion fits of these data 
were performed using a Gaussian broadening with a width of 40 fs. Our terminology is sueh that 
EDCs are energy distributions at fixed ky, not at fixed emission angle. 
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